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I. INTRODUCTION 

The most noted power law in economics is perhaps 
the Pareto law, first observed by Vilfredo Pareto p], 0, 
Q more than a century ago. It was found that in an 
economy the higher end of the distribution of wealth f(x) 
follows a power-law 



~ x 



(i) 



and a is an exponent (now known as the Pareto expo- 
nent) which Pareto estimated to be « 3/ 2. For the last 
hundred years the changes in value a ~ 3/ 2 in time and 
across the various capitalist economies seem to be small. 
This implies that there is inequality in the wealth dis- 
tribution, and only a few persons hold the majority of 
wealth. 

In 1931, Gibrat [4( suggested that while Pareto's law 
is valid only for the high wealth range, the middle wealth 
range is given by the probability density 
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where xq is a mean value and a 2 is the variance. The 
factor (3 = 1/ \Jl(j 2 is also know an as the Gibrat index, 
and a small Gibrat index corresponds to a uneven wealth 
distribution. 

An unequal wealth distribution is associated not only 
to these functions, but also to any other one which is 
not of the form a Dirac ^-function. The problem of the 
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appearance of inequalities seems therefore to be rather 
general and not necessarily related to a particular shape 
- e.g. the power-law form of the Pareto law - of the 
wealth distribution, despite distributions can vary from 
case to case assuming qualitatively different shapes. In 
fact wealth distribution has always been a prime concern 
of economics. Classical economists such as Adam Smith, 
Thomas Malthus and David Ricardo were mainly con- 
cerned with factor wealth distribution, that is, the distri- 
bution of wealth between the main factors of production, 
land, labor and capital. Modern economists have also 
addressed this issue, but have been more concerned with 
the distribution of wealth across individuals and house- 
holds. Important theoretical and policy concerns include 
the relationship between wealth inequality and economic 
growth. 

Wealth inequality metrics or wealth distribution met- 
rics are techniques used by economists to measure the 
distribution of wealth among the participants in a par- 
ticular economy, such as that of a specific country or of 
the world in general. These techniques are typically cate- 
gorized as either absolute measures or relative measures, 
and in the literature one of the most important debates 
is on the issue of measuring inequality. While one type 
deals with the objective sense of inequality, usually em- 
ploying some statistical measure of relative variation of 
wealth, the other type deals with some indices that try 
to measure inequality in terms of some normative notion 
of social welfare for a given total of wealth. 

Absolute criteria. Absolute measures define a mini- 
mum standard, and then calculate the number (or per- 
cent) of individuals below this threshold. These methods 
are most useful when determining the amount of poverty 
in a society. Examples include the poverty line, which is 
a measure of the level of wealth necessary to subsist in 
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a society. It varies from place to place and from time to 
time, depending on the cost of living and people's expec- 
tations. 

Relative criteria. Relative measures compare the 
wealth of one individual (or group) with the wealth 
of another individual (or group). These measures are 
most useful when analyzing the scope and distribution of 
wealth inequality. Examples include the Gini coefficient, 
which is a summary statistic used to quantify the ex- 
tent of wealth inequality depicted in a particular Lorenz 
curve. The Gini coefficient is a number between and 1, 
where corresponds with perfect equality (where every- 
one has the same wealth) and 1 corresponds with perfect 
inequality (where one person has all the wealth, and ev- 
eryone else has zero wealth). 

However, it has to be noted that "wealth" is here un- 
derstood differently respect to its common meaning: it 
represents the total amount of goods and services that a 
person receives, and thus there is not necessarily money 
or cash involved. Services like public health and ed- 
ucation are also counted in, and often expenditure or 
consumption (which is the same in an economic sense) 
is used to measure wealth. Thus, it is not clear how 
wealth should be defined. There is also the question that 
"Should the basic unit of measurement be households or 
individuals?" The Gini value for households is always 
lower than for individuals because of wealth pooling and 
intra- family transfers. The metrics will be biased either 
upward or downward depending on which unit of mea- 
surement is used. These and many other criticisms need 
to be addressed for the proper use of inequality measures 
in a well-explained and consistent way. 

In the attempt to answer some such basic questions 
and provide a foundation to the complex issues related 
to the appearance of wealth inequalities, various authors 
have independently formulated minimal models of wealth 
exchange [!, 0, 0, H, B G3 which, while being general 
enough to catch some universal features of economic ex- 
changes, are simple enough to be simulated numerically 
in detail and studied analytically. In these models a set 
of agents i = 1, . . . , M, representing individuals or com- 
panies whose state is defined by the respective wealth Xi : 
interact with each other from time to time by exchanging 
(a part of) their wealths. Such exchanges are defined by 
laws depending on the {xi} and also contain some ran- 
dom elements, as detailed below. A striking analogy was 
recognized - and actually motivated the introduction of 
some of these models - between the statistical mechanics 
of molecule collisions in a gas and these minimal models 
of economy, which are therefore referred to sometimes as 
kinetic wealth- exchange models. Such an analogy, rather 
than for its peculiarity, should be noticed since it signals 
a possible universal statical mechanism in action in the 
dynamical evolution of systems composed by single units, 
from gas composed of molecules colliding with each other 
exchanging their energy to economic societies in which 
single units interact by exchanging wealth. This analogy 
is being analysed here in more details than done previ- 



ously [Tl|, and represents the goal of the investigations 
presented here. 

We begin by recalling the main features of kinetic 
wealth-exchange models in Sec. |TT1 concentrating on an 
example of model with a fixed saving propensity X. 
In Sec. IIIII it is shown how the fact that for a sav- 
ing propensity A > one obtains an equilibrium 7- 
distribution j n (x) of order n, instead of the Boltzmann 
law ~ exp(— xj (x)) obtained for A = 0, actually strength- 
ens the kinetic analogy between economy systems and a 
gas in N(X) dimensions, with N(X) a known function 
of A [ll|, [l2[: through a general variational approach 
based on the Boltzmann entropy it is shown that the 
7-distribution j n (x) of order n is the equilibrium canon- 
ical distribution of a system with a quadratic Hamil- 
tonian H(qi, . . . ,qjsr) and N = 2n degrees of freedom 
(gi,..., qjy). The analogy is further discussed in Sec. HV] 
this time through a complementary microscopic approach 
based on the analysis of the dynamics of particle collisions 
in an TV-dimensional space. Through mechanical consid- 
erations only based on momentum and energy conserva- 
tion, it is shown that collision dynamics in N dimension 
can be recast in the form of the evolution laws of kinetic 
wealth-exchange models, both for A = and in the case 
with saving propensity A > 0, corresponding to a num- 
ber of effective dimensions N(X) > 2. Finally, in Sec. PVT 
conclusions are drawn. 



II. MAIN FEATURES OF KINETIC 
CLOSED-ECONOMY MODELS 

Simple social models of wealth exchange have been 
shown to well reproduce many features of real wealth 
distribution. For instance, the exponential law f{x) ~ 
exp(— f3x) observed at intermediate values of wealth is 
reproduced by a many-agent model system composed of 
M agents, who are assumed to exchange wealth in pairs 
at each iteration, according to the wealth-conserving evo- 
lution equations [13[ 

x ^ — e(^X"i I Xj ;) , 

x'j = e(xi +Xj). (3) 

Here e = 1 — e is a uniform random number in (0, 1), 
i and j are the labels of two agents chosen randomly at 
each iteration, and (xi,Xj) and (#-, x'j) represent the cor- 
responding wealths before and after a trade, respectively. 

More general versions of this model assign a (same) 
saving propensity A > to all agents, which represents 
the minimum fraction of wealth saved during a trade [5|, 
B GH Q HE • As an example, in the model of Ref. 
the evolution law is 

x\ = Xxi + e(l — X)(xi + Xj) , 

x f j = Xxj + e(l — X)(xi + Xj) . (4) 

It is to be noticed that while the total wealth is still 
conserved during a trade, x[+x f j — Xi+Xj, only a fraction 
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(1 — A) of the initial total wealth is reshuffled between the 
two agents during the trade. These models relax toward 
an equilibrium distribution well fitted by a 7-distribution 
7 n (x), as also noted by Angle [6], 



/r7(*) = 7n(0 

£ = fix, 



1 



T(n) 



e- 1 exp(-0, 



(5) 



where the scaling parameter is (3~ 1 = (x)/n and the pa- 
rameter 2n(A) = iV(A), as shown below, represents an 
effective dimension of the system and is explicitly given 



n(A) 



N(X) 



1 



3A _ 1 + 2A 
1-A ~~ 1-A 



(6) 



As the saving propensity A varies in A G [0, 1), the effec- 
tive dimension N(X) continuously assumes the values in 
the interval N G [1, 00). 

A comparison between results of numerical simula- 
tion of wealth exchange models and the 7-distribution 
is shown in Fig. HJ together with some real data for UK. 
The the zero-saving propensity model of Eqs. (|3]) is de- 
scribed by the exponential curve for A = 0. The main 
difference of the curves with A > respect to the expo- 
nential distribution is the appearance of a peak: while 
the average wealth is unchanged (the system is closed 
and one still has (x) = x tot /M, where M is the total 
number of agents) the number of agents with a a wealth 
close to the average value increase or, in other words, 
the wealth distribution becomes more fair for larger A's 
and therefore larger values of the effective dimension N. 
Eventually, as shown in Ref. f{x) — > S(x — (x)) for 
A —> 1 (N — > 00) thus becoming a perfectly fair distri- 
bution. As a consequence, measures of the inequality of 
the wealth distribution, such as the Gini coefficient, de- 
crease for increasing A and tend to zero for A — > 1 . In the 
rest of the paper we discuss the statistical mechanical in- 
terpretation of the equilibrium 7-distribution defined by 
Eqs. © and ©. 



III. MECHANICAL ANALOGY FROM THE 
BOLTZMANN ENTROPY 

It may seem that in going from the exponential wealth 
distribution obtained for A = to the 7-distribution cor- 
responding to a A > the link between wealth-exchange 
models and kinetic theory has been lost. In fact, 
the 7-distribution Jn/2( x ) represents but the Maxwell- 
Boltzmann equilibrium kinetic energy distribution for 
a gas in N dimensions, as shown in the Appendix of 
Ref. 0. Here a more general derivation of the 7- 
distribution Jn/2( x ) is presented, in which it is shown, 
solely using a variational approach based on the Boltz- 
mann entropy, that r )N/2( x ) is the most general canon- 
ical equilibrium distribution of an A/'-dimensional sys- 
tem with a Hamiltonian quadratic in the coordinates. 



Entropy-based variational approaches have been sug- 
gested (e.g. by Mimkes, Refs. [l8|, Oil) to be relevant 
in the description and understanding of economic pro- 
cesses. For instance, the exponential wealth distribu- 
tion observed in real data and obtained theoretically in 
the framework of the Dragulescu-Yakovenko model dis- 
cussed above, can also be derived through entropy con- 
siderations. Considering a discrete set of M economic 
subsystems which can be in one of J possible different 
states labeled with j = 1, 2, . . . , J and characterized by 
a wealth Xj, one can follow a line similar to the original 
Boltzmann's argumentation for the states of a mechani- 
cal system [18|, by defining the total entropy as 



W{rrij} = In 



Ml 



rai!ra2! . . . mjl 



(7) 



where rrij represents the occupation number of the j-th 
state; by variation of W{rrij} respect to the generic rrij, 
with the constraints of conservation of the total number 
of systems M = rrij and wealth x t ot = m j x ji one 
obtains the canonical equilibrium distribution 

rrij ~ exp(— (3xj) , (8) 

where f3 = 1/ (x) defines the temperature of the economic 
system. 

Here we repeat the same argumentation for an en- 
semble of systems described by a continuous distribution 
f(x), to show that this approach not only can reproduce 
the exponential distribution, but also the 7-distribution 
obtained in the framework of wealth-exchange models 
with a saving parameter A > 0, a natural effective di- 
mension N > 2 being associated to systems with A > 0. 

The representative system is assumed to have N 
degrees of freedom (gi, . . . , qjy) and a homogeneous 
quadratic Hamiltonian, that for convenience is written 
in the rescaled form 



q 2 n) = \q 2 , 



(9) 



where q = (#! + ••• + Qn) 1 ^ 2 ^ s the distance from 
the origin in the A^-dimensional g-space. As a me- 
chanical example, one can think of the N momentum 
Cartesian components (pi, . . . ,pn) an d the correspond- 
ing kinetic energy function K = (p\ + • • • + p%)/ 2m, 
where m is the particle mass; or of the Cartesian co- 
ordinates (#1, . . . , (Zjv) of an isotropic harmonic oscilla- 
tor with elastic constant n and potential energy U = 



K(q{ 



qj v )/2. It can be checked - e.g. using 



the Stirling approximation for the factorial function - 
that in the limit of large occupation numbers, the dis- 
crete version flTJ) of the Boltzmann distribution becomes 
proportional to — J dy f (y) In / (y) , where the continu- 
ous variable y replaces the discrete label j. For an N- 
dimensional system, the entropy will be proportional to 
-fdqi . - dq N f N {qi,...,q N )\n[f N {q 1 ,...,q N )}, where 
In(qi^ • • • 5 Qn) is the probability distribution in the 
A^-dimensional space. Then the constraints on the 
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FIG. 1: (Left panel) Comparison between numerical results (various symbols) of the equilibrium wealth distribution in the 
kinetic wealth- exchange models of closed economy and the corresponding gamma functions (continuous lines) with suitable 
parameters (see text for details). (Right panel) Distribution f(x) of individual weekly income or wealth x in UK for 1992, 1997 
and 2002; figure adapted from A. Chatterjee and B.K. Chakrabarti cond-mat/0709.1543v2 using data adapted from G. Willis 
and J. Mimkes, [cond-mat/0406694| 



conservation of the total number of systems reads 
J dqi . . . J dqN fN^Qi-, • • , Qn) — const and that on total 
wealth is j dq 1 . . . j dq N f N {q 1 , . . . , q N )X(q 2 ) x tot . In 
the end, using the Lagrange method, the equilibrium dis- 
tribution density can be derived by functional variation 
respect to /jvGzi? • • • •> Qn) of the effective functional 

Seff L/W] = J dqi... J dq N f N (q 1 , . . . , q N ) x 

x {ln[f N (q u . . . , q N )} + /i + f3X(q)} , (10) 

where \i and f3 are two Lagrange multipliers. Before car- 
rying out explicitly the variation it is convenient to sim- 
plify the problem exploiting the invariance of the Hamil- 
tonian, only depending on the variable q, under rotations 
in the g-space. The equilibrium probability density is ex- 
pected to depend solely on the distance q as well, i.e., 
/iv((Zi? n) = /jv(^)- This allows the reformulation of 
the variational principle as an effective one-dimensional 
problem in the variable q, by transforming Cartesian to 
polar coordinates, so that 



J dqi... J dq N (. . . ) -> j dq J 



dn (...), 



where the integration in Q is over the total solid angle in 
the TV-dimensional space. After integration of the N — 1 
angular variables, one is left with 



#eff [fl 



+co 

dqfi(q) 



In 



a{N)q 



N-l 



(11) 



where a(N) = / dQ = 2tt n / 2 / T(N/ 2) represents the 
surface of a hyper-sphere in N dimensions [20] with unit 
radius and we have expressed the probability density 
/iv((7i? • • • 7 Qn) — fN(q) m the A/'-dimensional space in 
terms of the reduced probability density fi(q) for the 
variable q, 



/ 



dnf N (q)= ( r(N)q N - 1 f N (q) 



Finally, moving from q to the energy variable x - 
X(q 2 ) = q 2 / 2, the corresponding probability density is 



/(*) 



dq(x) 
dx 



h(q)\ 



fi(q)\ 



q=q(x) 



q=q(x) 



2x 



where q(x) is obtained by inverting Eq. ©. In terms of 
the new variable x and distribution f(x), the effective 
functional (fTTj) becomes 



Sefflf] 



dx f(x) 



In 



a(N)x N /< 



-/i+/3x 



,(12) 



Variation of this functional, SS e ^[f]/6f(x) = 0, leads 
to the equilibrium 7-distribution ([5]) with dimensionless 
variable ^ = f3x and index n = N/ 2. As a simple ex- 
ample of application of this formula, one can obtain the 
Maxwell-Boltzmann probability density in three dimen- 
sions fs(K) for the kinetic energy K letting N = 3. In 
turn this suggests the interpretation of the parameter 
N (A) of wealth-exchange models as an effective dimen- 
sion of the system and of the Lagrange multiplier 
as the effective temperature, as it is in fact consistently 
recovered according to the equipartition theorem, 



2(s) 
N 



(13) 



IV. MECHANICAL ANALOGY FROM 
COLLISIONS IN N DIMENSIONS 

The deep analogy between kinetic wealth-exchange 
models of closed economy systems, where agents ex- 
change wealth at each trade, and kinetic gas models, in 
which energy exchanges take place at each particle colli- 
sions, can be further investigated and justified by study- 
ing the microscopic dynamics of interacting particles. In 
this section we make more rigorous an argumentation 
only mentioned in Ref. [ill ]. 

In one dimension, particles undergo head-on collisions, 
in which they can exchange the total amount of energy 
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they have, i.e. a fraction uj = 1 of it. Alternatively, 
one can say that the minimum fraction of energy that a 
particle saves in a collision is in this case A = 1 — uj = 0. 
In the framework of wealth-exchange models, this case 
corresponds to the model of Dragulescu and Yakovenko 
mentioned above [13], in which the total wealth of the 
two agents is reshuffled during a trade. 

In an arbitrary (larger) number of dimensions, how- 
ever, this does not take place, unless the two particles 
are travelling exactly along the same line in opposite 
verses. On average, only a fraction uj = (1 — A) < 1 
of the total energy will be lost or gained by a particle 
during a collision, that is most of the collisions will be 
practically characterized by an energy saving parameter 
A > 0. This corresponds to the model of Chakraborti 
and Chakrabarti [10[ , in which there is a fixed maximum 
fraction (1 — A) > of wealth which can be reshuffled. 

Consider a collision between two particles in an TV- 
dimensional space, with initial velocities represented 
by the vectors v (1) = . . . , V( 1)N ) and v (2) = 

( v (2)i ? • • 5 v (2)n)- For the sake of simplicity, the masses 
of the two particles are assumed to be equal to each other 
and will be set equal to 1, so that momentum conserva- 
tion implies that 



V(i) = v ( i) + Av , 



'(2) 



V( 2 ) - Av , 



(14) 



where and are the velocities after the collisions 
and Av is the momentum transferred. Conservation of 



'(2) 



'(1) 



energy implies that v^A + v 
using Eq. ([II]) , leads to 

Av 2 + (v (1) -v (2) )-Av = 0. 



v 2 2 ^ which, by 



and, furthermore, they can be considered as random vari- 
ables, if a hypothesis of molecular chaos is made con- 
cerning the random initial directions of the two particles 
entering the collision. 

However, the ?Vs are not uniformly distributed and 
their most probable value drastically depends on the 
space dimension: the greater the dimension TV, the more 
unlikely it becomes that the corresponding values (n) 
assume values close to 1 and the more probable that in- 
stead they assume a small value close to 1/N. This can 
be seen by computing their average - over the incom- 
ing directions of the two particles or, equivalent ly, on the 
orientation of the initial velocity v^) of one of the two 
particles and of the momentum transferred A v. In TV di- 
mensions, the Cartesian components of a generic velocity 
vector v = (v±, i>2, • • • , vn) are related to the correspond- 
ing hyper-spherical coordinates - the velocity modulus 
v and the (TV — 1) angular variables ipj - through the 
following relations, 

V\ = v cos (fi , 

V2 = v sin if i cos (f2 •> 



v N 



v sin cp i sin (f2 • • • cos <p n 



(19) 



(15) 



Using these transformations to express the initial velocity 
V(i) of the first particle and the momentum transferred 
Av in terms of their respective moduli and Av and 
angular variables {0^} and the expression (fTB|) for 

the cosine r\ becomes 

7*1 = cosai 

= COS 01 COS 0i + 

= [sin 0i cos 02] [sin 0\ cos 62} + 



Introducing the cosines n of the angles between the 
momentum transferred Av and the initial velocity v^) of 
the i-th particle (i = 1,2), 



COS C6j 



'(0 



V(i)Av 



(16) 



where = |v^)| and Av = |Av|, and using Eq. ([T5]h 
one obtains that the modulus of momentum transferred 



= [sin 0i sin 02 . . . cos 0at] [sin 9\ sin O2 . . . cos On] .(20) 

The average of the square cosine r\ is performed by first 
taking the square of ([20]) . integrating over the angular 
variables, considering that the TV-dimensional volume el- 
ement is given by 



N-l 



d N v 



~ idv n 



sin <ft j 



Av = -nv(i) J rr 2 v { 2) • 



(17) 



From this expression one can now compute explicitly the 
differences in particle energies X{ due to a collision, that 
are the quantities x\ — Xi = (v^X — v (i))/ 2. With the 
help of the relation (fT5|) one obtains 



x\ +r\x2~ r\ x\ , 
X2 -r\x2+ r\ x\ . 



(18) 



Comparison with Eqs. ([3]) for the kinetic model of Drag- 
ulescu and Yakovenko clearly shows their equivalence - 
consider that also here the r^'s are in the interval (0, 1) 



and finally by dividing by the total solid angle. In the 
integration only the squared terms survive, obtained from 
the square of ([20]) , since all the cross-terms are zero after 
integration - they contain at least one integral of a term 
of the form sin ip cos ip — sin(2y?)/2 which averages to 
zero. By direct integration over the angle 0^ and 0$, it 
can be shown that in TV dimensions 



(r 2 ) = (cos 2 ai) 



1 

TV 



This means that the center of mass of the distribution 
for 7*1, considered as a random variable due to the ran- 
dom initial directions of the two particles, shifts toward 
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smaller and smaller values as TV increases. The 1/N 
dependence of (rf) can be compared with the wealth- 
exchange model with A > 0. There a similar relation 
is found between the average value of the corresponding 
coefficients e(l — A) or e(l — A) in the evolution equa- 
tions (|4]) for the wealth exchange and the effective di- 
mensions N(X): since e is a uniform random number 
in (0,1), then (e) = 1/2 and from Eq. (|6|) one finds 
<e(l-A)> = (l-A)/2 = 6/(JV + 4). 

V. CONCLUSION 

The appearance of wealth inequalities in the minimal 
models and examples of closed economy systems consid- 
ered above appears to reflect a general statistical mecha- 
nism taking place for a wide class of stochastic exchange 
law - besides closed economy models - in which the state 
of the M units {i} is characterized by the values {xi} of 



a certain quantity x (e.g. wealth or energy) exchanged 
when units interact with each other. The mechanism 
involved seems to be quite general and leads to equi- 
librium distributions f(x) with a broad shape. In the 
special but important case of systems with a homoge- 
neous quadratic Hamiltonian - or equivalently with evo- 
lution laws linear in the quantities x\ - and N (effective) 
degrees of freedom, the canonical equilibrium distribu- 
tion is a 7-distribution j n (x) of order n = N/2. The 
corresponding distribution for the closed economy model 
with a fixed saving propensity A has the special prop- 
erty that it becomes a fair (Dirac-J) distribution when 
A — > 1 or N(X) — > oo. The possibility for single units to 
exchange only a fraction of their wealth during a trade - 
corresponding from a technical point of view to a wealth 
dynamics in a space with larger effective dimensions N 
- seems to be the key element which makes the wealth 
distribution less inequal. 



[1] V. Pareto, Cours d'economie politique, Rouge, Lausanne, 
1897. 

[2] V. Pareto, Cours d'economie politique, Reprinted as a 
volume of Oeuvres Completes, Droz, Geneva, 1971. 

[3] V. Pareto, Manual of political economy, Kellag, New 
York, 1971. 

[4] R. Gibrat, Les Inegalites Economiques, Sirey, 1931. 

[5] J. Angle, The surplus theory of social stratification and 
the size distribution of personal wealth, in: Proceed- 
ings of the American Social Statistical Association, Social 
Statistics Section, Alexandria, VA, 1983, p. 395. 

[6] J. Angle, The surplus theory of social stratification and 
the size distribution of personal wealth, Social Forces 65 

(1986) 293. 

|http://www .jstor.org] 

[7] E. Bennati, La simulazione statistica nell'analisi della 
distribuzione del reddito: modelli realistici e metodo di 
Monte Carlo, ETS Editrice, Pisa, 1988. 

[8] E. Bennati, Un metodo di simulazione statistica 
nell'analisi della distribuzione del reddito, Rivista Inter- 
nazionale di Scienze Economiche e Commerciali 35 (1988) 
735. 

[9] E. Bennati, II metodo Monte Carlo nell'analisi econom- 

ica, Rassegna di lavori dell'ISCO X (1993) 31. 
[10] A. Chakraborti, B. K. Chakrabarti, Statistical mechanics 

of money: How saving propensity affects its distribution, 

Eur. Phys. J. B 17 (2000) 167. 
[11] M. Patriarca, A. Chakraborti, K. Kaski, Statistical 

model with a standard gamma distribution, Phys. Rev. 

E 70 (2004) 016104. 
[12] M. Patriarca, A. Chakraborti, K. Kaski, Gibbs versus 



non-Gibbs distributions in money dynamics, Physica A 

340 (2004) 334. 
[13] A. Dragulescu, V. M. Yakovenko, Statistical mechanics 

of money, Eur. Phys. J. B 17 (2000) 723. 
[14] J. Angle, Deriving the size distribution of personal wealth 

from the rich get richer, the poor get poorer, J. Math. 

Sociol. 18 (1993) 27. 
[15] A. Chakraborti, Distributions of money in model markets 

of economy, Int. J. Mod. Phys. C 13 (2002) 1315. 
[16] J. Angle, The statistical signature of pervasive compe- 
tition on wage and salary incomes, J. Math. Sociol. 26 

(2002) 217. 

[17] M. Patriarca, A. Chakraborti, K. Kaski, G. Germano, 
Kinetic theory models for the distribution of wealth: 
Power law from overlap of exponentials, in: A. Chat- 
terjee, S.Yarlagadda, B. K. Chakrabarti (Eds.), Econo- 
physics of Wealth Distributions, Springer, 2005, p. 93. 
arxiv.org: physics/0504153 

[18] J. Mimkes, G. Willis, Lagrange principle of wealth 
distribution, in: A. Chatterjee, S.Yarlagadda, B. K. 
Chakrabarti (Eds.), Econophysics of Wealth Distribu- 
tions, Springer, 2005, p. 61. 

[19] J. Mimkes, Y. Aruka, Carnot process of wealth distribu- 
tion, in: A. Chatterjee, S.Yarlagadda, B. K. Chakrabarti 
(Eds.), Econophysics of Wealth Distributions, Springer, 
2005, p. 70. 

[20] E. W. Weisstein, Hypersphere. From Mathworld — A 

Wol fram web resource. 

://mathworld.wolfram.com/Hypersphere.html 



http 



